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Abstract. A new type of gauge fixing of the Coulomb gauge domain wall 
fermion system that reduces the fluctuation of the effective running coupling 
and the effective mass of arbitrary momentum direction including the region 
outside the cylinder cut region is proposed and tested in the 16^ x 32 x 16 
gauge configurations of RBC/UKQCD collaboration. 

The running coupling at the lowest momentum point does not show in- 
frared suppression and compatible with the experimental data extracted from 
the JLab collaboration. The source of the fluctuation of the effective mass 
near momentum p =0.6GeV region is expected to be due to the domain wall 
fermion zero modes. 

1 Introduction 

The infrared (IR) QCD is characterized by the color confinement and the chiral 
symmetry breaking. The features are studied via Dyson - Schwinger equation 
and/or via lattice simulation, but the interpretation of data does not necessar- 
ily agree with each other. We performed [Ij lattice simulations of the domain 
wall fermion which preserves the chiral symmetry in the zero - mass limit us- 
ing the full QCD gauge configurations of the RBC/UKQCD collaboration of 
16'^ X 32 X 16 lattice [2j and compared with results of staggered fermion of the 
MILC collaboration p9] . In these lattice simulations and in comparison with other 
works, we observed qualitative differences between quenched and unquenched 
simulations in the QCD effective running coupling, and in the Kugo-Ojima pa- 
rameter which is calculated numerially by the inverse of the Faddeev-Popov 
operator at the zero momentum [U O |6]. 

In O |8] we observed the running coupling of unquenched lattice simulation 
measured from the quark - gluon coupling and that measured from the ghost- 
gluon coupling are consistent, and they have enhancement as in comparison to 
the two loop perturbative QCD (pQCD) results in the region of the momentum 
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q around a few GeV. The running coupling in the IR hmits is suppressed in 
the quenched approximation [5]. Using the unquenched configurations in the 
Landau gauge, we observed also a suppression, however in the Coulomb gauge 
[8], we obtained the running coupling which shows no infrared suppression and 
consistent with the experimental data extracted by the JLab collaboration p]. 

In the Gribov-Zwanziger Lagrangian of the ghost-gluon system, there is a 
term that breaks the Becchi-Rouet-Stora-Tyutin(BRST) symmetry softly, 
which does not affect the renormalizability of the theory but affects the gluon 
propagator and the ghost propagator in the infrared [TO] . In the analytical calcu- 
lation, by adding a new mass term of the auxiliary field that is introduced in A^, 
the gluon propagator and the ghost propagator that are consistent with lattice 
simulation data are discovered. In quenched lattice simulation of relatively large 
lattices, the infrared suppression of the Landau gauge running coupling is con- 
firmed [H] , but it is due to the absence of the infrared singularity of the ghost 
propagator [121 E] • 

The difference of the quenched and unquenched lattice simulation in the IR 
region could be attributed to some topological effects. In the finite temperature 
Coulomb gauge QCD, van Baal pH] studied instantons (carolons) on X R. 

On the the Polyakov loop Pi{x) = —TrP exp{i / Ai{x -\- sei)ds} is defined, 

2 Jo 

and in the Aq = gauge by applying an anti-periodic gauge transformation, one 
could produce a non-trivial topology with the twist Pi{t = T) = —Piit = 0) 
and in the limit of T — > oo, instanton in the T"^ x R space was produced. The 
domain walls of RBC/UKQCD collaboration are separated by the distance of 
Lg in the 5 th dimension. In the limit of — > cxd the topology of the domain 
wall fermion and the gauge field may become diffeomorphic to x R. Although 
supersymmetry is not explicitly realized on the lattice, the fermionic zero modes 
could affect the topology of the bosonic gauge fields, as shown in the Lagrangian 
field theory [15]. 

In [T], we studied the domain wall fermion propagator with the cylinder cut, 
i.e. in the momentum region whose direction is close to that along the diagonal 
of the four dimensional lattice. Since the gauge configurations of RBC/UKQCD 
have the time axis twice as long as the space axis, we calculated the propagator of 
the cylinder cut which has the momentum components {p,p,p, 2p). In the region 
outside the cylinder cut, the fluctuation of the data of the effective running 
coupling and the effective mass were too large in the analysis using samples of 
the order of 100. We conjecture that the fluctuation is due to the fermionic zero 
mode and via a specific gauge fixing in the 5 dimensional space, based on the 
quaternion real condition, reduce the fluctuation. In this paper, we extend to 
some general momentum configurations. 

Organization of this paper is as follows. In sect. 2 we present the method 
of choosing the gauge in the 5 dimensional space and in sect. 3, we check the 
smallness of the deviation from the quaternion real condition, for the moment 
{Px,Py,Pz,Pt) = {p,p,p,0),ip, 0,0,0) and {p,p,p,2p) with p = pa = 1,2,3,4 (a 
is the lattice spacing). The validity of the Ansatz is checked Applications to the 
measurement of the effective running coupling and the effective mass are given 
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2 The gauge fixing in tlie 5tii dimension 

We specify the 5th dimension of the domain wah fermion by s = {0, 1, ■ ■ ■ , — 1} 
and define the fermionic part of the Lagrangian as [161 ITTl [T8] , 

Spill', 1p, U) = - ^ i>x,s{DF)^^s;y,s'1py,s', (1) 
x,s;y,s' 

where 

{DF)x,s;y,s' = "^s.s'^ly + Sx,yDj-si.. (2) 

The action Z?" contains the gauge field, and the interaction in the 5th dimension 
defined by D-^ does not contain the gauge field. 

Using the Dp defined in eql2]we make a hermitian operator Dh = j^R^Dp, 
where (i?5)ss' = 5s^Ls-i-s' is a reflection operator as 

/ -mp^PL 1->Pr 75(^" - 1) \ 

I^Pr 75(^11-1) I^Pl 



Dh 



I^Pr 75p"-l) I^Pl 
V 75(-D" - 1) 75-Pl -rufj^PR j 



where Pr/l = (1 ± 75)/2, and -1 in — 1) originates from D^^,. 
The quark sources are sitting on the domain walls as 

q(x)=PL^(x,^)^PR^{x,L,-\). (3) 

We take Pl1P'(x,s) oc e"^^)^, Pr^{x,s) oc e'^^^T"^)^, (s = 0, 1, • • • , - 1) with 
r = 0.842105, which corresponds to (r/Ls)/(l + r/Lg) = 0.05. 

The spinor basis that we adopt for the domain wall fermion is 75 diagonal 

i.e. 

a„, \ f I 



where 70 = ( j g ^ and o"j (i=l,2,3) are the Pauli matrices 

-CO' ''-(^o'). W 

The domain wall fermion propagator of the momentum p can be expressed 

as 

^ -i^ + M^p) -ii) + M{p) 

^p^ + M{p)M^{p) ^^ ^p^ + M^{p)M{py^^ ^ ' 

where using the plane wave with the coordinate of the 5th dimension s specified 
as x(P;*) s-iid the solution of the conjugate gradient method s), overlap 
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^{xiP,s)PL/R'^{p,s)) oc Bl/r and Tr(x(p, s)i;^-PL//j!^'(p, s)) oc Al/r, the mass 
is defined as 

^ Re[Bnip,Lj2)] 
Re[AR{p,Ls/2)] 

and 

Re[AL{p,Ls/2)y 

Here, {■p)M.[p) has a zero mode, while M.{p)M.\p) has no zero mode[T9]. 
In [T], we studied the domain wall fermion propagator in the Coulomb gauge 
adopting the cylinder cut, i.e. diagonal directions in the 4-d momentum space 
rix X riy X riz X rit spacial extension are selected. In this paper we consider other 
momenta. 

2.1 The U(l) real condition 

Since we adopt the representation of fermion spinors such that the 75 is diagonal, 
the fermion field ■0 can be gauge transformed in the 5th dimension as 

such that on the left and on the right domain wall, the propagator becomes 
approximately real by choosingflj 

be minimum (a- type), or 

be minimum (6- type), where we define 



rr(x(p,0)</'L(p,0)) 

\Tr{xip,0)MpM 
iOn ^ Tr{x{p,Ls - l)(f>R{p,Ls - 1)) 
\Tr{x{p,Ls-l)MP,Ls-l)\ 

Here, Tr is the trace in the color-spin space. 



(7) 



Tr{x{px,0/Ls - 1)(I)l/r{Px,0/Ls - l))ai + Tr{x{py,0/ Ls - l)(l)L/R{Py,0/ Ls - 1))(72 

+Tr{xiPz, - l)0^(p., 0/L, - l))a3 + Tr{xip,, - l)0^(p„ - 

= (a.°/^=-Vi + a,o/^^-V2 + a/Z^^-^as + at^/^^'-^lW detA^/^ (8) 

where 

The spinor of the fermion on the domain walls can be expressed in quaternion 
bases H which is obtained by adjoining an element j to the complex field C, i.e. 
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H = C + Cj ~ pO]. The equation in 4 dimensional space is converted in 
projective 3 dimensional space P3 and the antilinear map cr : P3 — P3 with 
cr^ = 1 defines a real structure on P3. 

The instanton solution proposed by Atiyah et al. [21] is produced in the 
quaternion bases and the transformation matrix which satisfies a specific con- 
dition which is cahed quaternion real condition was crucial for estabhsing the 
self-duality [2^ . We study whether the domain wall fermion on the left boundary 
can be transformed to that on the right boundary by this type of transformation 
matrix. 

In the lattice simulation of the domain wall fermion [Ij we sample wise mea- 
sure the propagator by using the transformation matrix from the left boundary 



which is given by the sample average of color diagonal components to the right 
boundary 



2.2 The quaternion real condition 

In the analysis of instantons in quaternion bases, Corrigan and Goddard |22j 
defined the transition function g(ij, vr) where vr is the complex two component 
spinor which satisfy. 



Here A is s non-zero complex constant, uj = pn where p is a quaternion and 
the transformation g{p7r, vr) between the domain walls should satisfy the reality 
condition [22]. In this case, a specific Ansatz can be introduced, which is discussed 
in sect. 3. 

©When Lj = p-K where p = p^ — ip-a, g{uj, it) is a quaternion can be expressed 

as 

g{p7r,7r) = h{p,Q)k{p,Cy^ 

where C = ^1 ^{P-, C) is regular in |(^| > 1 — e and /^(p, C,) is regular in |C| < 1 + e. 
In [22], the Ansatz 




was proposed as the transformation matrix. In our 5-dimesional domain wall 
fermion case, 7 = ^ — 1^ and ^, v contain the phase in the 5th direction ir/. 





g(Aa;, Avr) = (^(oj, vr), detg = l. 





(9) 



2/x = iuj2/T^2 -ir] = {px + iPy)C, + ipt - Pz - irj 



(10) 
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2v = icJi/vTi +irj= {pix - ipy)C + ipt + Pz + irj (11) 
The quaternion reality condition of the transformation matix 51(7, (^) gives 

OL.-i ^L.-i A ( f \^( ee-f f \(ao bo 

CL^.i dL^.i Mo C^e-f J V C-'eT M CO do 



(12) 



where / = /(7,-i). 

The function / and / taken in [22] is. 



In general cq and Ci^_i are polynomials of (, and they satisfy Ci^_iC^ = 
coC~^- We define. 

V' = c_iC"^ + ciC^ + 5 

where c_i = cLs-i and ci = co and 5 is a constant, which is defined later. 
The difference 

/Ar/fi'=^«^-i ( C'e-^ f \_( C'^-'' f \(ao bo 

' V cL.-i ^iL.-i^/^ v r^^y V r'e^yVco do 

(13) 

should be small, if the left wall and the right wall are correlated by the self- 
dual gauge transformation. With the Ansatz t\) = cqC"^ + CLs-iC^ + ^-^ we find 

C} = X — ^ — and the and the 5 are calculated from the simultaneous equation 

V CL.-l 

(HI. 



OL.-i/ + bL,-iC^e~^ = boC'e-^ + dof . . 

CL,^if + dL,^^C^e~^ = doC-^e\ ^ > 

The equation has two sets of solutions, which can be obtained numerically, using 
Mathematica[24J. 

Since 5(7, C)^o5~"'^(7) C) = ^Ls-i arid since the factor depends on the 
boundary conditions, there appears slight difference in the a-type minimization 
and the 6-type minimization. We found globally a-type gives smaller deviation 
than the 6-type, and we adopt the a-type in this paper. From the numerical 
practice, we observe that the deviation AL/R becomes small when the \5\ is 

large, i.e. when -0 ~ 5 ~ 4 is large. It means that e^"' ~ — ^ is a good 

^ do 
approximation in the case of small deviation. 



3 Numerical check of the deviation from the reaUty condition 

We calculated the quark propagator from 149 samples of the full QCD gauge 
configurations DWFoi produced by the RBC/UKQCD collaboration [2]. We first 
fixed the gauge to the Landau gauge and then to the minimum Coulomb gauge. 
We did not perform the remnant gauge fixing of Aq, since we do not want to affect 
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the topological structure of the ^o- The propagator calculation was done by using 
the conjugate gradient method [1]. We checked the deviation of the transforma- 
tion matrix using the phase rj of a-type U(l) real condition and parameter 6 and 
C of quaternion real condition, for the momentum {px,Py,Pz,Pt) = ip,P,P,'2p) 
i.e. cylinder cut, {p,p,p,0) i.e. no energy transfer and (p, 0,0,0) i.e. momentum 
transfer along the x-axis, using 49 samples of DWFqi configuration. 

We choose = ^J ^ to make {AL/R)2i consistent with and obtain e' 

and 5, which approximately satisfy the quaternion reality. The good indicator 
of the reality condition is {AL/ R)i2. In the tables, relatively large mean value 
and the same order of magnitude of the standard deviation means presence of 
exceptional samples. 

The results for the momenta in cylinder cut are shown in Table[TJ The devia- 
tion in the case of p = {p,p,p, 0) is shown in Table[2] and that for p = (p, 0, 0, 0) 
is shown in TableO 



p 


L/R 


iAL/R)u 


{AL/R)i2 


{AL/R)2i 


{AL/R)22 





L 


2.6(9.9) X 10"^ 


1.38(0.89) X 10-4 





1.14(2.18) xlO"'^ 




R 


3.4(10.4) xlO-*^ 


1.14(0.07) xlO""^ 





1.67(1.65) xlO"'^ 


1 


L 


2.4(4. 2)xl0-2 


0.77(1.56) 





0.17(0.48) 




R 


0.204(0.405) 


0.546(0.523) 





0.062(0.114) 


2 


L 


0.302(0.611) 


0.87(1.01) 





0.30(0.69) 




R 


0.98(4.61) 


1.59(5.81) 





0.96(4.93) 


3 


L 


0.72(1.62) 


1.26(1.94) 





0.65(1.18) 




R 


0.87(2.47) 


1.36(2.13) 





0.74(2.17) 


4 


L 


1.55(3.17) 


1.72(2.23) 





1.25(2.92) 




R 


0.51(1.20) 


0.91(1.28) 





0.54(1.61) 



Table 1. The deviation of the absolute value of the components of the transformation matrix 
for p — {p,p,p, 2p). Numbers in the brackets are the standard deviations. 



p 


L/R 


{AL/R)ii 


{AL/R)u 


{AL/R)2i 


{AL/R)22 


1 


L 


4.3(1.4)xl0"^ 


1.6(1.3) 





0.13(0.02) 




R 


5.9(l.l)xl0-3 


0.117(0.007) 





3.97(0.06) X 10-3 


2 


L 


0.62(1.47) 


1.44(1.89) 





0.56(1.08) 




R 


2.05(1.84) 


1.42(1.10) 





0.63(0.81) 


3 


L 


2.53(9.24) 


4.0(12.5) 





2.56(7.51) 




R 


0.52(0.52) 


0.73(0.61) 





0.58(0.47) 


4 


L 


1.52(2.90) 


1.80(2.82) 





1.64(2.35) 




R 


2.56(4.49) 


1.95(3.31) 





1.79(4.14) 



Table 2. The deviation of the absolute value of the components of the transformation matrix 
for p — {p,p,p, 0). Numbers in the brackets are the standard deviations. 



In the case of p = (4/o, 0, 0, 0) we observe an exceptional sample which has 
a large deviation, but we omitted this sample from the average. 
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p 


L/R 


{AL/R)u 


{AL/R)i2 


{AL/R)2i 


{AL/R)22 


1 


L 


0.214(0.314) 


7.18(4.68) 





1.11(0.94) 




R 


0.145(0.877) 


0.217(0.238) 





0.08(0.47) 


2 


L 


0.276(0.674) 


1.80(5.85) 





0.36(0.84) 




R 


4.67(4.77) 


1.49(1.17) 





0.42(0.70) 


3 


L 


0.303(0.441) 


12.1(7.0) 





2.2(1.3) 




R 


0.027(0.011) 


0.24(0.06) 





0.017(0.009) 


4 


L 


2.55(8.02) 


2.98(5.59) 





2.39(6.44) 




R 


3.84(8.07) 


1.84(5.26) 





2.18(5.99) 



Table 3. The deviation of the absolute value of the components of the transformation matrix 
for p — {p, 0, 0, 0). Numbers in the brackets are the standard deviations. 



4 The effective running coupling and the effective mass of the quark 

Corresponding to the sample-wise selection of larger absolute value of 6, i.e. 
large zero-mode component, we assign a parameter ind = 1 for larger and 2 for 
smaller, and multiply the phase (—1)*"'^ to the expectation value of the quark 
wave function that is used in the calculation of the effective mass of the quark. 
This phase factor can be absorbed in the phase e^^ introduced to approximately 
satisfy the U(l) real condition. 

We expect this phase ambiguity comes from the finite size of the length of the 
5th dimension, and this phase reduces the fluctuation dramatically. The mass 
function as a function of p(GeV) is shown in Figure 1. 



0.4 - 



3" 

^ 0.2 - 



0.0 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 



p(GeV) 

Figure 1. The effective mass of the domain wall fermion. niu ~ 0.01/a,ms = 0.04/a, without 
reality selection (upper dashed line) and with reality selection (lower dashed line). 



The fitted parameters of the mass function of the DWF and the staggered 
fermion are given in Table 2. In [25], the self-dual gauge configuration with 

= 1 had large fluctuation on configuration, and it was claimed that the 
Nahm transformation does not produce |Q| = 1 configuration in the system. 
The 5 dimensional space of the domainwall does not have periodicity in the 
5th dimension and the instanton like object that relates the fermion on the left 
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niud/O' ragja c yl(GeV) a 



DWFoi 


0.01 


0.04 


0.49 


0.76(2) 


1.25 


DWF02 


0.02 


0.04 


0.48 


0.80(3) 


1.25 


DWF03 


0.03 


0.04 


0.61 


0.66(2) 


1.25 


MILC/i 


0.006 


0.031 


0.45 


0.82(2) 


1.00 


MILC/2 


0.012 


0.031 


0.43 


0.89(2) 


1.00 



Table 4. The fitted parameters of mass function of DWF(RBC/UKQCD) and staggered 
fermion (MILC). 



wall and the right wall are not produced by the Nahm transformation on the 
fermionic zero mode. The self-dual configuration constructed a la ADHM and 
the Ansatz of Colligan and Goddard for establishing the quarternion reality was 
used in performing the gauge fixing of the domain wall fermion. The additional 
phase that is introduced following the sample dependent size of the zero-mode 
allows the instanton-like effects to become manifest. 

As shown in Figure 2, the QCD effective coupling of DWF in cylinder cut is 
slightly larger than that of the staggered fermion (dashed line) , but after applying 
the reality condition, it becomes slightly smaller. An enhancement in the region 
of g ~ a few GeV is expected to be due to condensates [271 EE]- Although it 
is necessary to study the continuum limit, the simulation of this size suggests a 
large part of the fermion configurations on the two walls are correlated by the 
self-dual gauge fields which plays the same role as instantons. 




-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
Log^o [q(GeV) ] 



Figure 2. Tlie effective coupling cvs.gi (?) of domain wall fermion and the gluon without reality 
selection (greenjand with reality correction(magenta) in cylinder cut region. Dark low momen- 
tum points are that with reality correction. Below 1.2GeV, red points are the experimental 
data of the JLab collaboration and the dark point is our result of momenta q — 0.7GeV/c. 
The dash-dotted line is the two loop pQCD result and the dashed line is the pQCD with 
condensates efffect derived from the fitting the running coupling [7] calculated by using the 
gauge configuration of fuU-QCD staggered fermion of the MILC collaboration [29] . 



As was shown in [T], the effective coupling in the IR region is experimentally 
investigated from the spin dependent Bjorken sum rule at large momenta and 
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Drell-Hearn-Gerasimov sum rule at low momenta [9]. In the literature there is 
also a Bethe-Salpeter analysis of the meson spectra [30]. Our data is consistent 
with the JLab data and suggest saturation of the coupling constant to a finite 
value and inconsistent with the meson data which suggest infrared suppression. 
The meson data at the lowest momentum point are based on transitions at high 
angular momentum states and rather ambiguous [31 j. 

5 Discussion and Conclusion 

The fermionic field on the domain wall at s = can be transformed to that 
on s = — 1 by a gauge transformation whose deviation from the quaternion 
reality AL/R is small. The quaternion real condition fixes parameters e^,C, and 
5 and we distinguish the two sets by the absolute value of 5, i.e. the parameter 
characterizing the zero-mode. One solution with large absolute value of the zero 
mode yields smaller AL/R. This selection reduces fluctuation of the effective 
mass of the quark and other physical quantities derived from IR - QCD domain 
wall fermion lattice simulation. 

It remains to check the finite size effects by extending the analysis to larger 
lattices [26]. It might be interesting to check how large the deviation from the 
reality condition AL/R is in the quenched configurations. 

Extension of the Nahm transformation which was done in the x R con- 
figuration [T3] to the X R and to study the instanton would be an interesting 
subject. Although there remains an explicit calculation of the field of the in- 
stantons, we have observed that there is an evidence of the self-dual gauge field 
contribution in the domain wall fermion propagator. 
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